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Planet formation: Statistics of spin rates and obliquities of 
extrasolar planets 
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ABSTRACT 

We develop a simple model of planetary formation, focusing our attention on 
those planets with masses less than IOMq and studying particularly the primordial 
spin parameters of planets resulting from the accretion of planetesimals and produced 
by the collisions between the embryos. As initial conditions, we adopt the oligarchic 
growth regime of protoplanets in a disc where several embryos are allowed to form. 
We take different initial planetary system parameters and for each initial condition, 
we consider an evolution of 2xl0 7 years of the system. We perform simulations for 
1000 different discs, and from their results we derive the statistical properties of the 
assembled planets. We have taken special attention to the planetary obliquities and 
rotation periods, such as the information obtained from the mass and semi major 
axis diagram, which reflects the process of planetary formation. The distribution of 
obliquities was found to be isotropic, which means that planets can rotate in direct or 
indirect sense, regardless of their mass. Our results regarding the primordial rotation 
periods show that they are dependent on the region where the embryo was formed and 
evolved. According to our results, most of the planets have rotation periods between 10 
and 10000 hours and there are also a large population of planets similar to terrestrial 
planets in the Solar System. 
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1 INTRODUCTION 

Following the first discovery of an extrasolar planet around 
51 Peg (jMavor fc Quelozlll995l ). the number of exoplanets 
known has risen to 429. Although most of them are gi- 
ant planets, the improvements in observational techniques 
have ensured that planets with masses less than 15Afe 
have started being detected with radia l velocity survey (e.g ., 
" 20061 ): lUdrv et al.l J2006h : iBonhls et all (|2007h ; 

and gravitational 



Lovis et al. 



Udrv et al.l 



: lUdrv et al 
?D; iMavor et all (1200 



microlensing survey l|Beaulier et al.ll2005l ). 

Although most of extrasolar planets so far discovered 
are giant planets, several statistical models for planetary 
growth presented in the last years suggest that a large 
number of small planets who fail to have enough mass to 
start the gas accretion onto the core exists dlda fc Linll2004l ; 
iMiguel fc Bruninill2009l : lMordasini et al.ll2009l ). and has still 
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not been able to be discovered (|Mordasini et all 120091 ) . At 
the time, several projects are in progress to detect terrestrial 
planets, we expect that they may And more Earth-size plan- 
ets in a close-future, but today, the sample is not enough 
and we also have to rely on what we know from our own So- 
lar System, and through computational models of planetary 
formation. 

This evidence supports the standard scenario, where 
terrestrial planets are formed through the next different 
stages: 1) agglomeration of dust particles through phys- 
ical collisions and setting in the protoplanetary disc, 2) 
planetesimal formation f r om grains in a thin m i dplan e 
|Goldreich fc Ward! 1 19731: IWeidenschilling fc CuziH 1 19931 ). 
3) runaway (e.g., [j Kokubo & Ida (Il996l) ) and oligarchic 
i|lda fc Makinolll993l : rKokubo fc Idalll998l ) accumulation of 
planetesimals to form protoplanets and 4) giant impact 
stage, where the embryos formed b y oligarchic grow th collide 
with one another to form planets (Wethcrill l 19851 ). 

The final stage of terrestrial planetary formation is the 
particular importance as it has a deep effect on the final 
characteristics of the planets: mass, orbital and spin param- 
eters. After this stage of planetary formation, the spin pa- 
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rameters of the planets change and evolve due mainly to 
tidal interactions with their satellite and host star. All of 
the terrestrial planets in our Solar System do not main- 
tain their primordial spin state and this is the reason why 
we unknown what primordial planetary spin would be ex- 
pected to find in a protoplanet. So questions as, what are 
the typical obliquity and rotation period that characterise 
the primordial planets? and how many collisions suffers a 
planetary embryo along its firsts years of formation? remain 
uncertain. 

A few works de aling with the study of p lanetary spins 
have been presented. iDones fc Tremaina (|l993l ) have exam- 
ined the accretion rate of spin angular momentum by a 
planet immersed in a differentially rotating disc of planetes- 
imals. They determined the mass and spin accreted by the 
embryos as a function of the velocity dispersion of the disc 
particles and the ratio of the planetary radius to the Hill 
radius. They found that if a protoplanet grows by accreting 
a large number of small planetesimals the spin angular mo- 
mentum of the planet will be determined by the called "or- 
dered component", but if a few giant impacts occur, most 
of the spin will be contribu ted by the "stochastic compo- 
nent" . lOhtsuki fc Ida! (|l998f ) have investigated the spin of 
a planet which accreted in a disc of planetesimals with non 
uniform spatial distribution. They results show that the or- 
dered component can dominate the final spin of the planet 
only if half of the size of the planet was acquire by the ac- 
cretion of small planetesimals and the size of the impactors 

is not too large. 

On the other hand lAgnor et al.l (1 19991 ) and IChambersl 

l|200ll ) have studied through N-body simulations the last 
stages of the terrestrial planet formation. They analysed the 
planetary obliquities as those found only considering the 
impacts between large embryos and have shown that this 
obliquities are expected to be represented by an isotropic 
distribution, result t hat was confirmed and generalized by 
iKokubo fc ldal(2007). who also considered an N-body code, 
but analysed a larger sample of embryos considering the 
standard disc model. 

Our principal aim is to make a statistical study of the 
primordial spin parameters of planets (obliquity and rota- 
tion period), resulting from the accretion of planetesimals 
and also due to the collisions between the emerging embryos. 
To this end we take different initial conditions, meaning dif- 
ferent discs, stars, initial number of embryos, and study the 
primordial planetary spins in different systems with the in- 
tention of obtain a better understanding of what we should 
expect to find in the Universe. We also analyse what are 
the consequences of planetary impacts in the mass and semi 
major axis diagram, considering embryos with masses less 
than 10M®. Our semi-analytical model takes as initial con- 
dition the oligarchic growth regime of protoplanets and al- 
lows them to migrate, fact that has a huge influence on the 
number of collisions suffered by an embryo. We adopt a per- 
fect accretion in collisio ns, supposition that was also consid- 
ered by other autho rs i|Agnor et al.lll999l ; [chambers! l200ll ; 
IKokubo fc Idall2007l) . but which says that the results should 
be interpreted cautiously. 

Each one of the 1000 systems considered, evolves for 
2xl0 7 years and we analyse the results statistically, finding 
an isotropic distribution of obliquities and where most of the 
planets rotate with a period between 10 and 10000 hours. 



We also found a large population of planets with the char- 
acteristics of terrestrial planets in the Solar System. 



2 MODEL AND BASIC EQUATIONS 

In this section we explain the model and basic equations con- 
sider in the work. As we take special attention to the plane- 
tary spin, the model adopted for the acquisition of angular 
momentum due to accretion and collisions between the em- 
bryos will be explain in detail. On the other hand the model 
for planetary growth and orbital evolutio n is essentially the 
same develo ped in our previous works iMiguel fc Bruninil 
(HH, 2009), which is a very simple model based on the oli- 
garchic growth regime and consider type I and 11 migration. 
For the sake of completeness, we will summarize it briefly 
below. 

2.1 Planetary growth 

We consider a protoplanetary nebula s tructure based on the 
minimum mass solar nebula (MMSN) (|Havashilll98ll) . where 
the surface density of solids at a distance a from the central 
star is 

Sd = 7/d7? lce (-^-) 7 gem' 2 (1) 

with rii ce a step-function which takes the value 1 inside the 
ice condensation radius and 4 outside it, expressing the effect 
of water ice formation. The snow line is located at ai Ce = 
2.7(jjj^J au from the central star of mass M*. 

On the other hand the volume density of gas is 

_ 3 

p gas (a, z) = p 9 ,o(a)e*W gem (2) 

where p gfi {a) = 1.4xl0" 9 / 9 (3^) ^S™" 3 - 

The parameters fd and f g state the solid and gas mass 
in the disc in terms of the MMSN model. We consider a 
large population of discs (1000 in each simulation), where 
we assume that f g follows a Gaussian distribution in terms 
of logiof g , centered at 0, with dispersion of 1 and fd is taken 
as fd = / 9 10 01 in order to consider more metallic discs. 

Both discs are not time-invariant. The gaseous disc 
change globally, decaying exponentially with a character- 
istic time-scale of Tdisc, which takes values between 10 6 — 
10 7 years in accordance to curre nt estimates of disc li f etimes 
around young solar type stars iBeckwith fc Sargentl (|l996l ) 
and the solid disc change locally, suffering the depletion of 
planetesimals produced by the effect of core's accretion. The 
disc of planetesimals also interacts with the nebular gas, this 
gas drag effect cause a radial motion of planetesimals before 
they become large enough to decouple f rom the disc gas 
l|Adachi et al.|[l976l ; iThommes et alj|2003h . we also consider 
this effect which was exp lained in detail in our previous work 
l|Miguel fc Brunini|[2009l '). 

The protoplanetary discs are extended between a» n — 
0.03442(|f|-) 2 au (|Vinkovidl2006l ) and 30au. The first initial 
core is located at a = a; n , the rest of the cores are separated 
10rH each other until the end of the disc is reached. Their 
initial masses are given by the minimu m mass necessary for 
starting the oligarch ic growth stage (|lda fc Makinol 1 19931 ; 
IKokubo fc Idall 19981 ). 
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with m the effective planetesimal mass. 

The solid accretion rate for a core in the oligarchic 
grow th regime, cons idering; the particle-in-a-box approxima- 
tion l|Safronovll 19691 ) is 



dM c 
dt 



10.53E d ft i?„ 1 + 



2GM t 



(4) 



where ft is the Kepler frequency, R p and Mt are the planet's 
radius and total mass (solid and gas) and a is the velocity 
dispersion which dep ends on the eccentricit y of the plan- 
etesimals in the disc. iThommes et al.l <|2003h obtain an ex- 
pression for the rms eccentricity of the planetesimals when 
gravitational perturbation of the protoplanets are balanced 
by the dissipation due to the gas drag, which is 



1.7 m 1 / 15 M\ lz pT 5 
6Vs C 1 / 5 p^o M 1/3 aV5 



(5) 



where b is the orbital separation between the embryos in Hill 
radius units, (6 = 10), Cd is a dimensionless drag coefficient 
which is ~ 1 and p m is the planetesimal bulk density. With 
this expression they found the next oligarchic-regime growth 
rate which includes the evolution of the planetesimal rms e 
and i, 



dM c 
dt 



^bTC^G^Mlplas^d 

— 1 J- 2 

Pm PM a 10 m 15 



Ml 



(6) 



where pM is the embryo bulk density, which is equal to the 
planetesimals density, then hereafter pu = Pm = p- 

The growth of the cores terminate when the solid sur- 
face density in their feeding zones is zero, which is caused by 
a combination of these factors: the embryos consume plan- 
etesimals on their feeding z ones, the density of planetesimals 
is dim inished by ejection l|Thommes et al.l 12003: Ida fc Lin! 
12004 ) and the planetesimal migration caused by the gas drag 
effect collaborate to empty this zone. 

Once the core became massive enough to retain a gas 
envelope, the effect of this atmospheric gas drag on the plan- 
etesimals increases the collision cross section of the pro- 
toplanet. This process was also taken into account in the 
model. 

When the core reaches the critical mass, the gas accre- 
tion process begins. In this work only those embryos with 
very few gas are considered, because the process of collisions 
between gas giant is poorly understood. For this reason we 
considered only those embryos with masses Mt < 10M®. 
Nevertheless I will explain the gas accretion model consid- 
ered for those protoplanets which attain the critical mass 
necessary to start the gas accretion process before reaching 
the 10M ffi . 

We assume that the critical mass necessary to start the 
gas accretion process is given by 



M, 



(- 



lQ-^M^yr- 1 ' 
This process occurs on a rate, 



dM g 
dt 



Mt 



(8) 



where M g is the mass of the surrounding envelope and r 9 is 
its characteristic growth time, 

mZ) yrs (9) 



this va lues were fitted fro m results obtained bvlFortier et al.l 
(2007) as is explained in lMiguel fc Bruninil (120081 ). 



2.2 Angular momentum transfer due to the 
accretion of planetesimals 

Our model also includes the acquisition of spin angular mo- 
mentum by the growing embryos due to the accretion of 
mass in the form of planetesimals. Mutual impacts between 
embryos contribute to the stochastic component of the an- 
gular momentum. On the other hand, accretion of a large 
number of small planetesimals produces an ordered spin an- 
gular momentum, which will be discuss in this section. 

In order to model the angular momentum accreted by 
the protoplanets d ue to the planetesim a l mas s accretion, we 
follow the work of iDones fc Tremainel (|l993h . Their model 
depends on two parameters: 

• The relevance of the velocity dispersion of the planetesi- 
mals in the planet's neighbourhood respect to the differential 
rotation of the planetesimal disc, and 

• the importance of the planet's gravity as compared to 
the self gravity of the disc. 

In the oligarchic growth regime, it is straightforward to 
demonstrate that the approp riate regime is that of high 
dispersion and strong gravity ()Dones fc Tremainelll993l ). In 
this case, if we analyse the contribution of the small plan- 
etesimals, we would found that the stochastic component is 
near one order of magnitude smaller that the ordered one. 
Therefore, we add to our model, only the ordered accretion 
of angular momentum due to the planetesimal accretion. 

According to the app ropriate three dimensional case of 
IDones fc Tremainel (|l993l ). the z component of the angular 
momentum L due only to the ordered component is given 
by 



Ll,ord — Mt ft 2 Rp ( yjA 



where 



A = 



R P ° 2 



RH = a{ M: 



1/3 



and the velocity dispersion is 

— — 2 o 2 2 

2 ^ ^ m 



(10) 



(11) 



(12) 



We assume that the RSM eccentricity of planetesimals 
in the disc is the equilibrium value found bv lThommes et all 
(|2003l ). which is given by equation Introducing equations 
Illlll2l and[5lin equation llOl we obtain the expression for the 
L z component due to the accretion of planetesimals, 



0.462 



M 



1/3 



m 2/15 p 3/5 



(13) 



Where all the units must be in cgs and the L z ,ord is in 
gem 2 Then at each time step, the z component of the 
angular momentum L changes by an amount 
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2/5 , ,1/6 ,,2/3 

AL z , ord ~ 6.5AxlO~ 7 p3a ^ I 1 AMt 



P 



11/15 a l/10 



(14) 



where AM t is the mass accreted in the form of plan- 
etesimals during the given time step. All the units are in cgs 
and the AIj. >orl j is in gem 2 s _1 . 

In order to obtain an estimate of which is the angular 
momentum acquired by a planet due only to the ordered 
component, we calculated the dependence of L z>or d with the 
embryo's mass for a protoplanet located at 1AU from a star 
like the Sun and we found, 



i(lAU) = 1.34a;10 4 



Mt 



5/3 



(15) 



Using equation 1131 we can obtain the rotation period 
reach by an embryo which only acquire angular momentum 
due to the accretion of planetesimals, which is, 

1/15 / v 1/10 / x -2/5 

/Vd-1501 



3g/cm 3 



a 
lau 



Pgas 



p gaa (lau) 



-1/6 



hours 



(16) 



for a planet located at 1AU and which orbits a star with 
IMq is Pord — 150 hours. 

Nevertheless this is not the only mechanism able to 
change the spin of the emerging embryos, the collisions be- 
tween the protoplanets have a huge importance in order to 
determine the final spin parameters. 

2.3 Collisions 

In the later stages of planetary formation, collisions repre- 
sent an important evolutionary process which plays a signifi- 
cant role in determining the final mass and spin state of the 
planets. These interactions are not fully understood, here 
we explain the model considered in the work, which is very 
simple but enables us to get some conclusions regarding the 
primordial obliquities and rotation periods of planets. 

When two protoplanets are too close to each other, mu- 
tual gravitational influence can pump up their eccentricities 
to values sufficient to ensure their orbits to cross. Once the 
protoplanets have perturbed one another into crossing or- 
bits, their subsequent orbital evolution is governed by close 
gravitational encounters and violent, highly inelastic colli- 
sions. 

Under the assumption of perfect accretion in collisions, 
we consider that a merger between protoplanets will occur 
if their orbital spacing, Aa, is less than 3.5 Hill radius. 

The magnitude of the relative velocity at which two 
bodies of total masses Mt,i and Mt,2 and radii _Ri and R2 
collide is 



Vcol = (Vrel + v\ ) - 



(17) 



where v re i is the relative velocity between the two bodies 
far form an encounter and v e is the scape velocity from the 
point of contact, given by 



V e = 2G 



Mt,i + M t . 



(18) 



\ Rj + R2 
The relative velocity between two embryos of orbital 
velocities v\ and V2 is 



Vrel 



Vx — V 2 



Considering that 02 = 01 + Aa, with Aa << ai and 
that the collisions are randomly oriented we obtained the 
following equat ion which sho ws the relative velocity between 
the embryos (| S afronovl 11969 ). 



Vrel 



Aa 



(19) 



with Q the orbital angular velocity (to our degree of ap- 
proximation it is equivalent to adopt a — a\ or 02, but we 
choose as a the semi major axis of the more massive planet). 
The distribution of velocities is isotropic, so the direction is 
chosen randomly with an isotropic probability distribution. 

We assume that in the beginning the embryos do not 
rotate but during its evolution they acquire spin angular 
momentum by the accretion of planetesimals (as seen in sec- 
tion [2]2} and by the collisions with other embryos. Here we 
analyse the total spin angular momentum of the resultant 
embryo acquired after a collision which is, 



L ir. 



(20) 



with L S pi n the sum of the spin of the target, L S pin,tar, and 
the spin of the impactor, L sp m,im and L co i is the spin angu- 
lar momentum delivered by the impactor during a collision 
where the impact point on the surface of the target is ran- 
domly calculated by assuming spherical embryos. 

Our assumption of perfect accretion occasionally allows 
particles to spin faster than break-up and destroy the em- 
bryo. This happen when the acceleration produced by the 
rotation are higher than gravity, which means, 

GM t 



Rlj 2 > 



R 2 



(21) 



with lj the rotation angular velocity. This condition leads to 
a critical value for the angular velocity, 



Ldcrit 



( GM t 
\ R 3 



(22) 



beyond which the embryo is gravitationally unbound. 



2.4 Orbital evolution 

When a protoplanet is embedded in a disc, their interaction 
are significant and lead to different regimes of planetary mi- 
gration regarding the embryo mass. When the protoplanet 
involved is a low-mass planet, the interaction can be calcu- 
lated using a linear theory which leads to a type I planetary 
migration, but when the planet reaches the mass necessary 
to open up a gap in its orbit, the disc response can no longer 
be treated as linear and it leads to the type II regime. The 
critical mass is derived from the condition 

r H > h (23) 

which is necessary for open a gap (|Lin fc Papaloizoulll993h . 
and where h is the disc scale of height. 

The mode l for type I and II migra tion is the same con- 
sidered before dMiguel fc Brunini 2009), which is essentially 
the same used by llda fc Linl (|2008r ) in their model, where 
the time-scales are given by 
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Figure 1. Histogram of the number of collisions suffered by each 
embryo at the end of the simulation. We see that a few collisions 
determine the primordial planetary spin. 
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with a = 10 -3 a dimensionless parameter which charac- 
terises the viscosity and the factor 



— is introduced for 

considering other important effects that might slow down 
the migration, without introducing a mayor degree of com- 
plexity to the model. 

We assume that both migration mechanisms stop when 
the core reaches the inner edge of the disc. 



3 RESULTS 

We investigate the statistical properties of primordial plan- 
etary spin resulting from the process of planetary formation 
through numerical simulations. In this section we show our 
main results. 



3.1 Some statistics of planets found 

We generate 1000 discs, for every system the mass of the 
star is taken random from values which follow a uniform log 
distribution in the range of 0.7 — IAMq, and the time-scale 
for the depletion of the disc of gas has a uniform distribution 
in log scale between 10 6 and 10 7 years. Each system evolves 
for 2a;10 7 years. 

We consider the formation of planetary systems which 
have suffered type I and II regimes of planetary migration, 
where the retardation constant for type I regime of migration 
is taken as C mg u = 0.1. 

Figure [1] shows an histogram of the number of collisions 
suffered by each embryo over the 2a;10 7 years. We note that 
most of the planets suffer less than 5 impacts during its for- 
mation, which means that in most of the cases primordial 
spins of planets are randomly determined by a very few im- 
pacts suffered during accretion. On the other hand we also 
found some planets which have more impacts. This is due 
to the migration of the embryos which makes some of them 
move rapidly towards the star and suffer more collisions than 
the most external ones, which have very few embryos to col- 
lide with. 

We also analyse the distribution of planetary primordial 
obliquities. Figure [2] shows this distribution where we see 
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Figure 2. Histogram of primordial obliquities found, which cor- 
respond to an isotropic distribution. 



that the obliquity distribution corresponds to an isotropic 
distribution of the spin vector, given by 

1 



p(e) = — sin(e) 



(26) 



this r e sult confirms th e ear lier findings o f lAgnor et al.l 
|l999fl ; IChambersi (|200ll ) and lKokubo fc Idal (|2007r i. which 
were obtained using N-body simulations, and is due to the 
fact that during this stage of planetary formation, the scale 
of height of the disc is much larger than the size of the 
embryos, so collisions can occur in any direction. For this 
reason, the result is indeed independent of the orbital evolu- 
tion of the embryo: the isotropic distribution is maintained 
if we consider planetary migration or if we do not, the only 
difference is the amount of collisions suffer by the embryos. 

Since the ordered component has obliquity of or 180 
degrees, we would expect that this will have consequences 
in the obliquities distribution, but the effect of the ordered 
accretion is important in those embryos which do not suf- 
fer any collision during their growth. In these cases, the 
angular momentum acquired by the accretion of planetes- 
imals determines their final spin state. On the other hand 
for those embryos that suffer great collisions during their for- 
mation, their obliquities are determined by the momentum 
acquired during the impacts because this stochastic com- 
ponent is very strong and dominates the final state of the 
embryo. According to our results, most of the embryos suf- 
fer collisions during their formation and for this reason we 
do not found significant changes in the obliquity statistics 
made. 

We also note that there are a large amount of embryos 
that did not collide with any other. These are tiny embryos 
who grew in a low mass disc and are located near to the inner 
edge of the disc, as a consequence, they were not eaten by a 
larger embryo migrating towards the star. 

We show in figure [3] the obliquity of the planets plot 
against their mass. This plot shows that we can found, 
with equal probability, terrestrial planets with obliquities 
between and 180 degrees, which means planets who ro- 
tate in a direct or indirect sense, independently of its mass. 
This result tells us that the primordial spins of planets are 
not those commonly observed in the terrestrial planets in 
our own Solar System, whose current spin axes are more or 
less perpendicular to their orbital planes (except for Venus) . 
However, the spin axes of the terrestrial planets are not pri- 
mordial, so this does not necessarily indicate a problem in 
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Figure 3. Obliquity of the surviving planets plotted against their 
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Figure 4. Distribution of rotation periods of all the planets found 
in our simulations. 

the model consider here. Other studies s uch as planet-hos t 
star or plane t -satel lite's tidal interaction iGoldreichl (l966); 
lAtobe fc Ida! (|2007n . among others, must be taken into ac- 
count for explaining the present obliquities of the terrestrial 
planets. 

3.1.1 The study of the rotation periods 

The rotation periods of planets were calculated assuming 
that the protoplanets were spheres of uniform density and 
the distribution found is shown in hgure [H where we note 
that most of the planets reach rotation periods larger than 
~ Whours but there are also a large amount of planets with 
periods between 0.1 and Whours. The planets who reach 
spin periods less than 0.1 hours are really rare, because at 
that rotation periods they use to have a angular velocities 
larger than the critical one. 

In order to understand the rotation periods distribu- 
tion, we study the angular momentum of the planets formed. 
Figure [5] shows the angular momentum as a function of 
the mass of the embryos, where the solid line represent the 
L C rit beyond which the planets are gravitationally unbound, 
which can be deduced from equation 1221 The dotted line 
shows the angular momentum due to the ordered compo- 
nent that acquire a planet located at 1AU from the Sun, 
which was deduced from equation [13] and shown in equation 
1151 We note that those embryos who did not suffer any col- 
lision and acquire angular momentum only by the accretion 
of planetesimals, should have L near the dotted line, while 
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Figure 5. Mass and angular momentum of the planets found 
in our simulations. The solid line represents the L cr u beyond 
which the embryos are gravitationally unbound. The dotted line 
represent the value for the ordered angular momentum reached 
by the cores located at 1AU from the sun (see equation 1156 . 



those who experienced the change of momentum due to one 
or more impacts, could reach a larger angular momentum 
but always below the stability limit. 

We know that the rotation period is inversely propor- 
tional to the angular momentum, so the higher the angular 
momentum, the shorter the period, as a consequence there 
can not be planets with small periods, and that is why we 
have an absence of planets with periods less than ~ 0.5 hours 
in the figure^] 

We also study the evolution of the rotation periods in 
three different simulation times: at 1000, 10 5 and 2a;10 7 
years, which is the final simulation time. 

The rotation periods of embryos as a function of its mass 
in the three different times is plot in figure [6] In the first 



figure (6(a) I, when 1000 years have passed, we found only 
embryos with rotation periods until 100 hours, and small 
masses, which is probably due to the short time that has 
passed, the embryos did not have much time to grow. We also 
observe a small population of more massive embryos with 
shorter periods, some of them with masses of up to 7M^. 
Since very few time has passed, they are probably embryos 
located in the interior region of the disc. This region is rich 
in solids and this favors the rapid growth of the embryos 
cores', which makes them the firsts to suffer a large amount 
of collisions and hence increase their spin angular velocities, 
leaving them to the brink of instability. 



As time passes (figure 6(b) I we note that the embryos 
acquire larger periods, and the amount of embryos with 
small periods decreases. Finally, at the end of the simula- 



tion (figure 6(c) I, we note a well-marked difference between 
the few planets with periods less than ~ 1 hours and the rest 
of the population. These are very rare planets. As seen in 
equation 1221 those embryos with small rotation periods ro- 
tate rapidly, so their spin angular velocities are high enough 
to overcome the critical rotation angular velocity for rota- 
tion instability. As a consequence we find a small amount of 
planets with this periods, only a very few percent survive, 
and the surviving ones have mainly small masses. 

On the other hand we observe that most of the planets 
have reached rotation periods of up to ~ 10000 hours. These 
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Figure 6. The rotation period is plot against the embryos' mass 
for different times. Figure [6(a) j presents the results at 1000 years, 
in figure [6(b) | has pass 10 years and the last figure, [6(c) | shows 
the distribution at the end of the simulation (2xl0 7 years). 



are probably the embryos that only acquire their angular 
momentum by the accretion of panetesimals. 

In figure [7] the rotation period is plot as a function of 
the embryo's semi major axis, where those planets with the 
largest rotation periods probably acquired them mainly by 
the accretion of planetesimals, while those with the shorter 
periods need one or more impacts for having that spin. 

We can also compare our results with those observed in 
the terrestrial planets in our own solar System. In the case 
of Mercury and Venus, whose current rotation periods are 
58.64 and 243.01 days respectively, their spin rates have un- 
dergone great changes since their formation and could not 
be considered as primordial. The close proximity of these 



Figure 7. The figure shows the rotation period and semi major 
axis of all the surviving planets at the end of the simulation. 



planets to the Sun produces a tidal dissipation that has 
slowed down the spin rates, then the primordial values of 
the rotation period must have been much lower than those 
currently observed, for this reason we can not compare with 
these planets. 

The case of Mars and the Earth is different, because 
they are located far from the Sun and that is why Solar 
tides have not altered their spins appreciably. The rotation 
speed of Mars can be considered as primordial, because its 
satellites are so small that have not influenced appreciably in 
the spin rates' evolution. In the case of the Earth, while the 
magnitude of the Earth-Moon system' s angular momen- 
tum has been approximately constant since its formation, 
the Earth' spin has been slowed down by lunar, since tidal 
interactions have transferred angular momentum from the 
Earth to the Moon. 

While the current (and primordial) rotation period of 
Mars is ~ 24.5 hours, the one that should have had if all its 
mass were obtained only by the accretion of planetesimals 
(no collisions involved), could be deduced from equation 1161 
and is 247.9 hours. In the case of the Earth, the maximum 
value reached by the rotation period since it was formed is 
the current one, and the one obtained only by the ordered 
component is 150 hours. According to this results Mars and 
the Earth did not acquire their rotation periods only by the 
accretion of planetesimals, but during one or more impacts 
during its formation. 

We also note that in the population of planets shown 
in figure [7] we found a large sample of planets with the 
characteristics of the Terrestrial Planets. 



3.1.2 Mass and Semi major axis distribution 

In our previous works iMiguel fc BruninU (|2008l . 120091 ). we 
have studied the changes in the ma ss and sem i majo r axis 
diagram due to different factors. As llda fc Lir] (|2004l ) have 
shown, this diagram shows the process of planetary forma- 
tion, where different regimes of planetary growth were found 
depending on the material available in the region on the disc 
where it was formed. 

As equation [5] shows, the cores' planetesimal accretion 
rate depend on the region where the embryo is located and 
the solids available, which are larger at the smallest (a < 1 
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au) semi major axis (equation [T|. For this reason, we found 
a rapid cores' growth in the inner regions of the disc, where 
a < 1 au, and the lowest solid accretion rates are found in 
the outer regions of the disc, where the embryos take longer 
to grow. 

As seen in previous section, the region in the disc where 
a planet grows, has also a strong influence in the rotation 
period reached by the planet, because a smaller semi ma- 
jor axis ensures that the planets have more solids available 
and they will accrete more angular momentum. Besides, if 
they are in a region with a large density of embryos, they 
would have a large probability of collisions, that will change 
their spin too. Figures |8(a)|8(b)| and |8(c)| show the mass 
and semi major axis distribution in the three times studied: 
1000, 10 5 and 2xW 7 years respectively. As seen in the first 
figure, in the begging the embryos who grow faster are those 
initially located really close to the star, in the inner part of 
the disc and as seen in figure [7] between these embryos we 
also found those with the smallest rotation periods. On the 
other hand, those embryos located in the intermediate and 
outer part of the disc remain almost with their initial mass. 



As time passes (figure 8(b) |, we note that the embryos who 
grow faster have migrated closer to the star or fragmented 
and disappear, while the embryos located in the intermedi- 
ate region began to grow. The results obtained at the end of 
the simulation, represented in figure 8(c) show that those 
embryos who were located initially in the intermediate re- 
gion of the disc, have migrated to the star (or fragmented 
and disappear) and those in the outer regions of the disc, 
have now starting to grow reaching rotation periods until 
- 10000 hours. 

In our previous works we have focused our attention to 
those planets with masses larger than 1M®. Here we study 
planets with masses less than IOMq and found another im- 
portant effect which changes the distribution of mass and 
semi major axis of terrestrial planets. As seen in before, em- 
bryos with small periods, can not resist the acceleration due 
to rotation and destroyed themselves. This fact will change 
slightly the mass and semi major axis distribution obtained 
in our previous works. With the aim of comparing with our 
previous results, we plot in figure O the mass and semi ma- 
jor axis distribution found when fragmentation by collisions 
was not taken into account. 

Comparing figure 8(c) with figure O we observe fewer 
planets considering the fragmentation by collisions that 
those found in the other case. So this is an important ef- 
fect that must be considered when working with terrestrial 
planets. 
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Figure 8. Mass and semi major axis distribution found in differ- 
ent evolution times. Figure [8 (a) | shows the results at 1000 years, 
in figure |8(b)| has passed 10 s years and in figure |8(c)| we plot 
the results at the end of the simulation. In all the simulations we 
consider the fragmentation of embryos by collisions. 



3.2 Those who did not survive 



As we have shown in section 12.31 there is a stability limit 
beyond which the planets are not able to remain united 
and disarmed. Here we show some statistics regarding those 
planets who could not survive. 

Figure [TOl shows the mass and rotation period of these 
"broken" embryos as they were when reached the critical 
rotation angular velocity. Figure |10(a)| presents the results 
of those embryos fragmented before the firsts 1000 years, 
who are mostly those with periods less than 0.1 hours, in 
figure [Tb(b)| we see the results at 10 5 years and finally the 
last figure (figure 10(c) I shows the total embryos who do not 



survive. As seen in the figures all the broken embryos have 
rotation periods less than ~ 2 hours, which is approximately 
the critical period. We also observe those embryos with small 
spin periods are the firsts fragmented, then the embryos with 
periods near to 1 hours, and finally those of ~ 2 hours 
exceed the limit of stability. 



4 SUMMARY AND CONCLUSIONS 

In the process of planetary formation protoplanets collide 
with one another to form planets. We have investigated the 
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Figure 9. Mass and semi major axis distribution found without 
considering the fragmentation of planets by collisions. 

final assemblage of terrestrial planets from protoplanets us- 
ing a simple model which consider the oligarchic growth 
regime of protoplanets as initial condition in a disc where 
several embryos are allowed to form. As explained in our pre- 
vious work the formation of several cores simultaneously in 
the disc has a strong influence on the dynamic of the plan- 
etesimal disc, which influences directly the growth of the 
embryos' cores and the final assemblage of planets found. 

In our model we also have included the interaction be- 
tween the protoplanets and the disc, which leads to a plan- 
etary migration. When a embryo is migrating towards the 
central star it could perturbate the cores placed in its path, 
causing the accretion of the core in most cases, this collisions 
affect the spin state of the embryos. 

As collision among giant planets are poorly understood 
we have focused our attention on planets with masses less 
than IOMq where a very simple model for planetary im- 
pacts has been considered. We suppose that when two em- 
bryos are a distance less than 3.5 Rh the merger between 
both protoplanets occurs, which leads to the union of two 
embryos to form a single body. This perfect accretion model 
produces spin rates that are too high and when the accel- 
eration produced by the rotation is greater than those of 
gravity the body overcome the critical spin angular veloc- 
ity for rotational instability and is fragmented. This simple 
model allows us to obtain some interesting results regarding 
the final properties of terrestrial planets. 

We also have considered the acquisition of angular mo- 
mentum due to accretion of planetesimals. The accretion of a 
large amount of planetsimals produces an ordered spin that 
adds angular momentum to that acquired during collisions, 
so the final spin of the planets is a result of this two effects. 

In order to analyse the statistical properties of the 
assembled planets we take different initial planetary sys- 
tem parameters, considering 1000 different discs, where each 
planetary system evolves 2xl0 7 years. 

As in our previous works we have analysed the infor- 
mation provided by the mass and semi major axis diagram, 
which reflects the process of planetary formation. We ob- 
serve fewer planets with masses less than 1M® considering 
the fragmentation by collisions that those found in the with- 
out this effect. This means that the effect of fragmentation 
by collision has a strong influence on the final population 
of terrestrial planets formed and should be considered when 
these planets are involved. 
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Figure 10. Mass and rotation period of embryos fragmented. 
Figure |10(a)| shows those embryos fragmented before the firsts 
1000 years, in figure |10(b)| we see the results at 10 5 years of 
simulation and figure [10(c) | shows the total embryos who do not 



We also have studied the effects produced by the col- 
lisions between the embryos, where we find that most of 
the planets suffer less than 5 impacts during its formation, 
which means that in most of the cases primordial spins of 
planets are randomly determined by a very few impacts suf- 
fered during accretion. 

We also take special attention to final spin state, which 
means planetary obliquities and rotation periods, where 
we found that the distribution of obliquities of final plan- 
ets is well expressed by an isotropic distribution, result 
that confirms thos e obtained previously by other authors 
lAgnor et all l| 1999ft : iKokubo fc Idal (|2007ft and is indepen- 
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dent on the planetary mass. This fact is in marked contrast 
to the terrestrial planets in our own Solar System, whose 
current spin axes are more or less perpendicular to their 
orbital planes (except for Venus). However, the spin axis 
of the terrestrial planets strongly depends on the gravita- 
tional perturbations from the other planets of the Solar Sys- 
tem that create a large chaotic zone for their obliquities. So 
all of the terrestrial planets could have experienced large, 
chaotic variations in obliquity in their history, and this is 
why their obliquities can not be considered as primordial 
|Laskar fc Robutell Il993l ). So the fact that the terrestrial 
planets in our Solar System present obliquities ~ 0° does 
not necessarily indicate a problem with the model consid- 
ered here. Other studies such as body and atmospheric tides 
and core-mantle friction among others, must be taken into 
account for explaining the present obliquities of the terres- 
trial planets. 

Regarding the findings on the rotation period, we found 
that the primordial rotation periods of terrestrial planets 
are dependent on the semi major axis, which means on the 
region where the embryos were formed and evolved. 

On the one hand we note a very small population of 
planets with small rotation periods (less than ~ 0.5 hours), 
which are very rare planets, because at that rotation periods 
the spin angular velocities are high enough to overcome the 
critical rotation angular velocity for rotation instability. 

On the other hand there are a large population of em- 
bryos with rotation periods until ~ 10000 hours. These 
planets with large rotation periods probably acquired them 
mainly by the accretion of planetesimals, while those with 
shorter periods need one or more impacts for acquire that 
spin. 

Another important result is that we have found a large 
population of planets with the characteristics of the Ter- 
restrial Planets, and our results suggest that they did not 
acquire their rotation period only by the accretion of plan- 
etesimals, but during one or more impacts during their for- 
mation. 
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